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ABSTRACT 

The accuracy of position measurements on stellar targets with the future Space Interferometry 
Mission (SIM) will be limited not only by photon noise and by the properties of the instrument 
(design, stability, etc.) and the overall measurement program (observing strategy, reduction 
methods, etc.), but also by the presence of other "confusing" stars in the field of view (FOV). 
We use a simple "phasor" model as an aid to understanding the main effects of this "confusion 
bias" in single observations with SIM. This analytic model has been implemented numerically in 
a computer code and applied to a selection of typical SIM target fields drawn from some of the 
Key Projects already accepted for the Mission. We expect that less than 1% of all SIM targets 
will be vulnerable to confusion bias; we show that for the present SIM design, confusion may be a 
concern if the surface density of field stars exceeds 0.4 star/arcsec^. We have developed a software 
tool as an aid to ascertaining the possible presence of confusion bias in single observations of any 
arbitrary field. Some a priori knowledge of the locations and spectral energy distributions of the 
few brightest stars in the FOV is helpful in establishing the possible presence of confusion bias, 
but the information is in general not likely to be available with sufficient accuracy to permit its 
removal. We discuss several ways of reducing the likelihood of confusion bias in crowded fields. 
Finally, several limitations of the present semi-analytic approach are reviewed, and their effects 
on the present results are estimated. The simple model presented here provides a good physical 
understanding of how confusion arises in a single SIM observation, and has sufficient precision to 
establish the likelihood of a bias in most cases. We close this paper with a list of suggestions for 
future work on this subject. 



Subject headings: Data Analysis and Techniques 

1. Introduction 

The Space Interferometry Mission (SI1V0) is be- 
ing designed by NASA's Jet Propulsion Labora- 
tory to provide a facility-class instrument for mea- 
suring the positions and proper motions of stars at 
optical wavelengths with micro-arc-second (/las) 
precision. This represents an improvement by sev- 



^also currently called SIM-PlanetQuest 



eral orders of magnitude over the precision of all 
existing astrometric instruments. For faint sources 
(V ^ 15 mag.), SIM will also be more than a fac- 
tor 10 better than any other future planned space 
mission, and will therefore uniquely permit many 
new classes of problems to be addressed. Such 
problems include: the direct measurement (for the 
first time) of the masses of earth-like planets in 
orbit around nearby stars; determining the dis- 
tances to stars by direct triangulation over the 
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whole Galaxy and out to the Magellanic Clouds; 
measuring the transverse motions of galaxies in 
the Local Group; and, establishing the shape of 
the dark matter distribution in the Galaxy. A de- 
scription of the instrument and its current science 
program is available at the JPL/SIM web sited 
The mission is now at the end of the detailed de- 
sign phase. After more than 15 years of develop- 
ment, all the major technical questions have been 
answered. New devices have been invented in or- 
der to provide metrology internal to the spacecraft 
at a level of a few tens of picometers, a fraction of 
the inter-atomic distance in a molecule of oxygen. 
The next major step is to begin construction of 
the instrument. 

SIM will be the second optical interferometer 
in space devoted to astrometry, following the Fine 
Guidance Sensors (FGS) on the Hubble Space 
Telescope (HST). However, SIM is a Michelson in- 
terferometer using separated collectors, quite dif- 
ferent from the filled-aperture white- light shearing 
interferometer design of the FGS. SIM includes 
three long-baseline interferometers housed on a 
common truss, each formed by two « 0.3 m aper- 
tures which compress their light beams and guide 
them through delay lines to beam combiners. Dur- 
ing operation, two interferometers are used for 
precision guiding of the spacecraft, and the third 
views the target of interest. Data are then accu- 
mulated by tracking the target until enough pho- 
tons have been recorded to achieve the particular 
science goal. 

The precision with which astrometry can be 
done on a specific stellar target with SIM will 
be limited by photon noise, the design and sta- 
bility of the instrument, and by the data calibra- 
tion and processing. We have some control over 
the instrument properties, which depend on the 
specific choices made when implementing the de- 
sign in hardware. Also, the operation and cal- 
ibration of the instrument can be optimized so 
as to maximize the achievable precision. How- 
ever, there is another source of error which may 
be present, and which is largely out of our con- 
trol; it does not reduce the ultimate precision with 
which a given single measurement can be made 



^See \http://planetquest.jpl.nasa.gov/3IM/simAndex.cfrr^ 
for descriptions of the current set of Key Projects. Almost 
half of the available 5- year Mission time is still unallocated. 



with SIM, but it may reduce the final accuracy 
of that measurement. This source of error arises 
because of the presence of other stars in the SIM 
field of view (FOV) which can "confuse" any sin- 
gle observation made on the target star. The light 
from these extraneous stars perturbs the measure- 
ment, so that the measured target position can 
differ from the true position. The difference is a 
bias which can reach a level of many times the 
single-measurement precision estimated from the 
instrument parameters alone. It is this single- 
measurement confusion bias which concerns us in 
this paper. 

It is important to emphasize that the final ac- 
curacy of the astrometric parameters (position, 
parallax, and proper motion) determined by SIM 
for any given target star will be a result of carry- 
ing out a complex program of several single mea- 
surements on that target, plus repeated measure- 
ments on many other stars for the determination 
of calibration and baseline orientation p arameters 
(|Bodenll2004lMilman fc Turvshevl26o3l ). Since we 
are concerned only with possible bias in a single 
measurement^ the details of the entire observing 
program are not directly relevant here; however, 
it also means that we can not quantify the con- 
sequences of confusion bias to the final accuracy 
with which the "end-of-mission" astrometric pa- 
rameters can be obtained on any specific target. 
It is clear that the effects of single-measurement 
confusion bias will generally diminish as more ob- 
servations are combined. But this also means that 
projects which involve only a few observations of 
a target (e.g. "targets of opportunity" , single par- 
allax measurements on nearby bright stars, etc.) 
may have a greater susceptibility to confusion bias. 

Specific aspects of confusion in astrometric 
measurements with SIM have bee n considered by 
severa l authors in the recent past. iDalal fc Griest 
(j200ll ) showed that the characteristic response of 
SIM's fixed-baseline interferometer as a function 
of wavenumber and delay can be used to refine a 
model of the distribution of confusing stars in the 
FOV. This model can then be used to correct the 
measured position of the target of interest, and in 
many cases where the level of confusion is not too 
great the astrometric accuracy can approach the 
measurement precision. Dalai and Griest success- 
fully applied their method to models of confused 
fields in the LMC in which « 16 faint stars are 
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scattered over the FOV around the « 19 magni- 
tude target star. Photon noise is included in these 
models. These authors then go on to consider 
the additional complication if one of the stars in 
the FOV changes brightness owing to a micro- 
lensing event, and show that an extension of their 
fitting algorithm to include the precision photom- 
etry provided by SIM's detectors permits even this 
apparently-intractable case to be handled almost 
as well. However, their method fails when the 
angular separation between any pair of sources 
in the FOV (as projected on the interferometer 
baseline) corresponds to a delay difference of < 2 
coherence lengths for the full bandpass of the de- 
tection system. This is a projected separation of 
25 milli-arcseconds (mas). Indeed, this is a general 
limit for SIM observations. Our approach is some- 
what simpler than that of Dalai and Griest, and 
yields some improvement in the minimum angular 
separation which can be measured, but the basic 
limitation can not be overcome. We will compare 
our approach to theirs in more detail in a future 
paper. 



Raiagopal et al.l (|200ll ) also considered a num- 



ber of specific cases of confusion on SIM astrom- 
etry. These authors introduced a graphical anal- 
ogy using phasors as an aid to understanding how 
errors in the target position arise from confusing 
sources in the FOV. Typical target fields were con- 
structed on a simulated image with grid spacing 
of 5 mas, and the amplitude and phase of the 
fringe which would be measured with SIM for a 
given wavelength on that image field was com- 
puted with a Fourier transform. Diffraction ef- 
fects at the edges of the (presumed w 1" square) 
SIM FOV were included in constructing the model 
image, and vector averaging of the individual (nar- 
row) SIM wavelength channels was used to simu- 
late the 1-dimensional apodization of the fringes 
over the FOV caused by the decreasing coher- 
ence of the frin g es as the bandwidth increases. 
Raiagopal et al. (l200lh were particularly inter- 
ested in modeling the effects of mispointing of the 
FOV in subsequent visits to the same target field 
when target proper motions were being measured; 
in that case the actual distribution of field stars 
changes as some disappear from one side of the 
FOV and others appear at the other side. They 
included photon noise, and also addressed the is- 
sue of how the size of the FOV defined by the field 



stop influenced the level of confusion. Their source 
field models were constructed to simulate SIM ob- 
servations of the position and proper motion of 
target stars in M31, the LMC, and the Galactic 
bulge. They concluded that the confusion-induced 
errors in position can often be significant (several 
times the measurement precision) for faint target 
stars but the proper motion errors are likely to 
be small. The errors are, as expected, smaller for 
wider measurement bands. 

In o ne other confusion-related study. ^Takeuch i et al 
(|2005l ) considered the effect of circumstellar disks 
on the measurement of stellar wobble during ob- 
servations aimed at detecting extra-solar planets. 
Their models showed that neither the motion of 
the disk mass center nor the contamination by 
disk light is a serious threat to detecting planets 
around pre-main sequence stars; the basic rea- 
son for the insensitivity of these observations to 
confusion from circumstellar disks is that interfer- 
ometers tend to resolve such an extended source, 
reducing its influence on the astrometry of the 
parent star. 

The studies summarized above have shown that 
confusion poses limits to the accuracy of any single 
SIM measurement, and have therefore succeeded 
in raising our awareness of this problem. How- 
ever, the detailed design of SIM has changed since 
those studies were done, and many of the changes 
will affect on the modeling results. The size of the 
collector and its central obscuration, the entrance 
aperture field stop defining the geometrical FOV, 
the transmission efficiency of the optics, the fringe 
disperser design (which defines the bandwidth and 
central wavelength of the spectral channels), and 
the QE and spectral response of the detector are 
now all much better defined. Furthermore, pre- 
vious studies have focused on specific science pro- 
grams which were already suspected to be pushing 
the capabilities of the instrument; they have not 
provided us with any general "tools" for under- 
standing and recognizing confusion, or for dealing 
with it. Previous studies have also often taken a 
statistical approach which is less suitable for an- 
swering direct questions about specific fields, such 
as: is a SIM observation of this particular target 
embedded in that particular field of stars likely to 
be confused? And, can the observation be done 
in such a way so as to reduce the confusion bias? 
What a priori information about the target and 



3 



the field would help? And, if the observation has 
already been taken, can we identify the effects of 
confusion in the data? These questions have pro- 
vided the motivation for the work described in this 
paper. 

This paper is organized as follows. In the 
next section, we summarize the basic Michel- 
son interferometer response as it applies to SIM. 
We then recall t he ph asor model introduced by 



Raiagopal et al.l (|2001r ) and elaborate upon it as 



a tool for understanding the behavior of confusion 
in SIM astrometry. Using this analytic model, to- 
gether with updated knowledge of SIM's instru- 
ment parameters, we have constructed a simula- 
tion code for evaluating the likelihood of confusion 
bias in any specific field; details are presented in 
Section [3l In Section IH we present single measure- 
ment confusion bias as a function of magnitude dif- 
ference and projected separation of an additional 
star present within the SIM FOV. In Scction[5l we 
apply this semi-analytic model to a number of tar- 
get fields drawn from the Key Projects which have 
already been chosen for the initial SIM science pro- 
gram. From this experience we then consider how 
the single-measurement confusion bias might be 
reduced through the addition of other information. 
The most useful additions appear to be knowledge 
of the approximate locations and spectral energy 
distributions (SEDs) of the target and of the most 
troublesome confusing stars in the SIM FOV. Fi- 
nally, in Section [21 we describe the limitations of 
our current approach. These limitations are pri- 
marily related to the simplified model of the focal 
plane of SIM which we have used here. In par- 
ticular, in this paper, we have not modeled the 
detailed mechanism by which the spectral disper- 
sion is implemented H nor have we considered the 
pixellation of the focal plane by the CCD detector. 
We have explored these points with the aid of a 
more elaborate model that includes these effects, 
and find that the biases estimated using this more 
elaborate model differ only by small amounts from 
those provided by the approach described hereQ 
We have therefore chosen to present the main is- 
sues relevant to SIM confusion with a minimum of 
complication, and leave the discussion of the more 



thin prism disperser turns SIM into an objective prism 
spectrograph on the CCD detector. 
^However, consideration of this more sophisticated model 
does suggest additional ways to reduce any confusion bias. 



elaborate instrument model to a future paper. 

Binary stars will be an important class of tar- 
gets for SIM, and are the topic of one of the major 
Key Projects. In these cases, the two stars are in a 
bound orbit and are physically close to each other. 
Typical binaries to be studied with SIM will have 
separations from about a few mas to 1000 mas and 
orbital periods from a few days to several years. 
Stars in crowded fields can occasionally mimic the 
effects of binaries if their projected separations 
become small for particular baseline orientations, 
but the effects of confusion from such "apparent" 
binaries can be reduced (or even eliminated) by 
rotating the interferometer baseline and repeating 
the observation. However, for "real" binaries, ro- 
tation of the baseline is an integral part of the mea- 
surement process. The goal of the binary obser- 
vation is to obtain the characteristics of the orbit, 
and this is done by measuring the positions of the 
components for a number of baseline orientations. 
These targets are sufficiently specialized that we 
have removed them from the list of crowded-field 
problems treated in this paper. Such observations 
treat binaries as "signal", whereas here we treat 
them as "noise" . A discussion of astrometry on 
binary stars with SIM is planned for a future pub- 
lication. 

2. Astrometry v^fith SIM 

Historically, the angular separation between 
two stars on the sky has been determined using 
telescopes of various sizes mounted on mechani- 
cal structures including large metal segments with 
regular rulings on their edges, or a combination of 
such devices with measurements of transit times. 
An astrometric interferometer converts the prob- 
lem from a measurement of angular increments to 
a measurement of a distance interval which can 
be determined very precisely using laser metrol- 
ogy inside the instrument. Thanks to the absence 
of atmospheric instabilities in the space environ- 
ment where SIM will operate, an additional level 
of precision worth several orders of magnitude can 
be added by directly measuring the phase of the 
fringe pattern using target photons. The main fea- 
tures of how these measurements are made with 
SIM will now be described. 
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2.1. The SIM interferometer 



The response pattern P of a Michelson in- 
terferometer as appUed to astrometry and imag- 
ing in astronomy can be written, under quasi- 
monochromatic conditions, as foUows: 

P{S,9,X)= (1) 
Po{l + Asin{f[5-B9])}. 

where, Po is the total Ught collected by the two 
collectors, A is the fringe amplitude, S is the in- 
ternal path delay, B is the baseline, 9 is the angle 
on the sky, and A is the mean wavelength. See Ap- 
pendix [X] for the derivation of this equation. For 
convenience we often use the mean wavenumber 
k = 1/A. 

The response of the interferometer as a func- 
tion of 9 projected onto the two-dimensional sky 
is sketched in Figure [1] Figure [2] shows the same 
interferometer response, now as a function of the 
total optical path difference A = 6 — B9 and A. 
This Figure is also known as the "fan diagram". 
In this diagram, stellar targets appear at specific 
values of total delay along the x axis corre- 
sponding to their actual positions 6i in the sky 
and the setting of the internal delay d, and so can 
be represented as a set of vertical lines. An ex- 
ample is shown in Figure [2] for a single target at 
A = +lfj,. Moving vertically in this diagram from 
shorter to longer wavelengths at a given delay A, 
we see that the fringe phase changes linearly if the 
delay is close to the location of an isolated tar- 
get, and will jump by ±27r radians if the delay is 
large. If the delay setting corresponds to the ac- 
tual location of the target, the fringe phase will be 
a constant; this is the position at which the total 
optical path delay difference between the two arms 
of the interferometer is zero. It is this particular 
property of the fringes which enables us to identify 
the position of the target precisely. 

2.2. How SIM v^forks 

Consider a single target star located at 6i near 
the center of the FOV as shown in Figure [3] (ignore 
the other stars located at 62 and 63 for the time 
being) . It is now convenient to consider the inter- 
nal delay to be made up of two parts, a "coarse" 
part and a "fine" part, such that 6 = 6c + Sf. The 
measurement begins by setting Sf = and adjust- 
ing Sc until fringes are found. The bandwidth is 




Fig. 1. — The response of the interferometer as a 
function of the angle on the sky, for a fixed inter- 
nal delay S and wavelength A. A dense pattern 
of 1-D fringes showing the response of the instru- 
ment as a function of position is here projected 
onto a 2-D FOV on the sky. The interferome- 
ter baseline orientation is indicated by the thick 
double-headed horizontal arrow. The amplitude of 
the fringes is apodized by a 1-D "coherence FOV" 
function which depends chiefiy on the bandwidth 
of the channel, as described in Appendix [Xj 
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Fig. 2. — The "fan diagram" ; a contour plot of the 
Michelson interferometer response as a function of 
total path delay A = (5 — B9 and wavelength. The 
fringes appear to 'fan out' at longer wavelengths, 
hence the name. The possible location of a target 
star is shown at A = 
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then progressively widened (by further binning the 
channels) and the coarse delay further adjusted so 
as to maximize the fringe contrast. A limit will 
reached when the delay difference Sc—Si is approx- 
imately within the coherence length appropriate 
for the maximum available instrument bandwidth. 
Defining 60 = B9q as this "final" setting ofSc, the 
goal of the subsequent steps in the measurement 
process is to measure Si— So = B{9i —9o) by scan- 
ning the fine delay in a series of small steps to a 
maximum of ~ ±1A around zero, and recording 
the fringes which appear in each of n narrow-band 
channels. The response function of equation[T]now 
reads: 



P{Sf,4>n,kn) = 

Po{n){l + An sm[2Trk^Sf + , 



(2) 



.)]} 



where </)„ — 27TknB{9o — 9i). A fit is then done 
to the data in each channel, yielding the fringe 
parameters Po{n), An, and (/)„. 

A plot of the resulting "phase spectrum" (/)„ 
vs. the mean wavenumber /c„ is shown in Figured] 
This plot will show jumps of the fringe phase if the 
delay offset 61 — So is larger than a typical fringe 
period. Taking the phase at one of the channels 
as a reference, the phases of the other channels 
can be "unwrapped" by adding or subtracting 2n 
when a discontinuity is encountered. A plot of the 
"unwrapped" phases will be a straight line par- 
allel to the X axis which we call the "delay spec- 
trum" , as shown in Figure [5] The final delay offset 
(and hence the angular position 6*1) can be com- 
puted from a simple average of the data in this 
Figure, although several other possibilities exist 
(jMilman k Basingej[20o3 ). 

We have glossed over a large number of calibra- 
tion issues, including the precise measurement of 
delay increments along the internal delay line, the 
precise length of the baseline, the calibration of the 
relation So ^ 9o, etc. In addition, the measure- 
ment we have obtained is an angle with respect 
to the direction of the interferometer baseline, but 
that direction is not yet known in any external ref- 
erence frame such as equatorial coordinates. Pre- 
cise determination of the baseline vectors used for 
all SIM observations is clearly an extremely impor- 
tant part of the whole astrometric program. These 
calibration and baseline determination problems 
are a major part of the SIM project. However, 
confusion bias acts as a perturbation on a single 




Fig. 3. — A potentially-confused FOV consisting 
of a target star at 9i and two field stars at 6*2 and 
93. 9o corresponds to the current setting of the in- 
ternal "coarse" delay (5c- Scanning the "fine" delay 
Sf and fringe fitting will lead to the identification 
of a position which can be different from the de- 
lay offset corresponding to the "true" location of 
the target. Thus the presence of field stars can 
introduce a bias in the measurement of the target 
star's position. 
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Fig. 4. — The phase spectrum, showing the fringe 
phase (j) vs. wavenumber. Note the "phase-wrap" 
jumps. 
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Fig. 5. — Delay spectrum estimated from the un- 
wrapped fringe phases of Figure |4] as described in 
the text. The average delay is 1309.02 nm, and 
corresponds to an offset of 30000.58 /las. The 
"glitches" on this figure are a result of simplifi- 
cations we have made in our numerical model for 
the dispersion of the fringes with wavelength. 



2.4. Visualizing confusion 

One can visualize the nature of confusion with 
the help of Argand diagrams from the theory of 
complex variables (e.g. jPhilips 1961) since, after 
subtraction of the total power term, the sine- wave 
fringe pattern of Equation [T] produced by an in- 
terferometer observing a single un-confused target 
star maps one-to-one with a vector in the complex 
plane having modulus equal to the fringe ampli- 
tude and argument equal to the fringe phase. Such 
vecto rs are called phasors{e.g. lGaskilllll978l : iHechtl 
20021 ). and the results of adding many fringes to- 
gether from many stars in a crowded field are eas- 
ily understood in terms of a vector sum of these 
phasors as shown in Figure [6l a dapted from Fig- 
ure 1 of iRaiagopal et al.l (|20011 ). The bright tar- 



SIM observation, so it is possible to discuss the 
origin and nature of confusion bias in single mea- 
surements without a detailed consideration of the 
entire astrometric program of SIM. 

2.3. Confusion bias in SIM measurements 

The operation of SIM as described in the pre- 
vious section makes the assumption that the FOV 
contains only a single target star, and this is ex- 
pected to be the case for more than 99% of the 
fields to be observed. However, the real sky may 
occasionally contain other stars in the foreground, 
background, or even associated with, the target 
star, as shown in the sketch of Figure [3l Since the 
fringe patterns of all the stars in the FOV are the 
same in the same wavelength channel, the result- 
ing fringe pattern for several stars together will 
mimic that of a single star, but with a different 
amplitude and phase. If the bandwidth is wide, 
and the (baseline-projected) positions of the con- 
fusing stars are well away from the target star, 
their influence will be attenuated. Nevertheless, it 
is clear that the presence of such extraneous "field 
stars" in the FOV can lead to a bias in the mea- 
sured fringe phase, and therefore to a bias in the 
measured position of the target. The final mea- 
surement precision may not be affected, but the 
final accuracy of the target star position will be 
reduced. 



Y 




Fig. 6. — Phasor diagram at a specific wavenum- 
ber k showing the vector addition of the complex 
visibilities of a bright target star, located at the 
origin of the coordinate system, with the visibil- 
ities of 10 fainter stars, located elsewhere in the 
FOV. 

get is assumed to be located at the origin of the 
(local) coordinate system, and therefore has zero 
fringe phase; it is represented by a vector along 
the real axis with some amplitude. The field stars 
have non-zero fringe phases; they are represented 
by shorter phasors with relative magnitudes and 
directions depending on their locations in the field 
with respect to the target. The field stars form a 
"noise cloud" at the tip of the strong phasor rep- 
resenting the astrometric target of interest, and 
their combined effect is to alter the final fringe 
amplitude and phase. The shift (p of the fringe 
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phase from the true position of the target phasor 
is the bias owing to confusion. Since the phases of 
the confusing sources change with observing wave- 
length, and at a rate which depends on their (pro- 
jected) displacement from the nominal field center, 
the small phasors all rotate at different rates and 
directions as the observing wavelength changes, so 
that the net confusion bias (f> = (f>{X) of the tar- 
get of interest changes with the wavelength of the 
observations. It is now obvious that the effects 
of confusion will depend on the spectral energy 
distributions (SEDs) of the target and of all the 
confusing sources. 

3. Simulating SIM's response 

In order to proceed with a simulation we need 
to know the expected structure of the distribution 
of light in (and surrounding) SIM's FOV, and the 
response of the instrument to this distribution. 

3.1. Source models 

In general, we assume that source distribu- 
tions can be approximated as a set of unresolved 
sources. We ignore the presence of extended emis- 
sion, which will in any case contribute little to the 
net fringe amplitude (but may increase the pho- 
ton noise on the detector) . Each of these sources 
is represented by a delta function with some am- 
plitude and some 2-D location in the FOV. The 
spatial distribution is either idealized from exist- 
ing images of the target area, e.g. images from the 
cameras of the Hubble Space Telescope (HST), or 
defined ad hoc with an educated guess. A model 
consists of a set of x, y positions within (and just 
outside of) the FOV, along with the brightnesses 
and the spectral energy distributions (SEDs) of 
the sources at each position. If the actual SEDs 
are not available, model SEDs are either approxi- 
mated or obtained from spectral libraries □ 



® Prom the web site 

^ttp://garv,et.stsci.edu/3TIS/stis^models.html#models' 



when spectral types are known and from the file 
kurucz^orig.fits available at the web site 

\http://dae45Aaa.csic.es:8080/'^jmaiz/software/chorizos/chorizos.ht^^ j-QllQ}lj)'lll 
when gravity, effective temperature and metallicity are 
known. 



3.2. Instrument models 

We provide here a high-level description of the 
SIM instrument model we have adopted in order 
to determine the total power, fringe amplitude and 
fringe phase of a given source distribution. The 
numerical details and relevant figures are collected 
in Appendix IB. 21 

3.2.1. Incident photon fluxes 

The number of photons detected from each star 
within the FOV is estimated as follows: 



Channel response 

We estimate the mean wavelength and bandpass of 
all 80 channels of SIM from the knowledge of dis- 
persion (wavelength vs. position at the focal plane) 
and detector pixel size. We assume that each pixel 
corresponds to one channel, and that the total en- 
ergy of all the stars within the FOV is collected in 
a single pixel. This is an idealization, but we defer 
details of the dispersion and the pixellation in the 
focal plane to a later paper where a more complete 
instrument model is presented. As we shall discuss 
further in fJ51 these additional complications have 
only small effects on our results, and come at the 
expense of considerable additional complexity. 

Apodization through the field-stop 

For each star in the FOV, the effects of diffraction 
through the entrance aperture and apodization by 
the field stop are modeled as the convolution of 
a circular stop with the (wavelength-dependent) 
diffraction pattern of SIM's entrance aperture in- 
cluding the effects of the central obscuration (cre- 
ated to accommodate the internal laser metrology 
system) . This gives the aperture transmission fac- 
tor (a fractional number as defined here) which, 
when multiplied by the total power, gives the frac- 
tion of the total light from the star that falls within 
the field stop. See Appendix IB. 21 for the details. 



The optical train of SIM involves many reflections, 
each with some loss which depends on the wave- 
length, and ultimately the light is imaged on a de- 
tector which itself has some wavelength response. 
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These various reflectivities have been multiphed 
with the net CCD detector quantum efficiency in 
order to obtain an overall wavelength-dependent 
"throughput" of the optical system. We have 
fitted this function to a polynomial for conve- 
nience in numerical computations; the coefficients 
are listed in Appendix lB.il The SEDs, originally 
available for a mag 10 star, are modified to ac- 
count for the actual magnitude of the stars and 
are multiplied by twice the collecting area of a 
single aperture, the aperture transmission factor, 
and the throughput, integrated over the channel 
bandpass, and expressed as photons/sec/channel. 

3.2.2. Fringe amplitude for each star 

The visibility amplitudes of the confusing stars 
are attenuated by the bandpass function (cf. Equa- 
tion |A6| according to their angular distances from 
the position of the target measured parallel to 
the baseline, and the wavenumber and coherence 
length for the specific channel. 

3.2.3. Fringe phase for each star 

The visibility phases (/)/ of the field stars de- 
pend on their projected angular distances 0f and 
are given by — 2TTkB sm{p cos{PA ~ = 
27rkB6f, where p is the radial distance of the field 
star measured in radians, PA is the position angle 
of the field star, and tp is the position angle of the 
baseline orientation. These quantities are defined 
in a Cartesian coordinate system with origin at the 
center of the FOV, and with the X and Y axes ori- 
ented along the directions of Right Ascension and 
Declination. 

3.3. Resultant fringe models 

The resultant fringe is the sum of fringes pro- 
duced by the target and each of the field stars. In 
the notation introduced for equation [21 this total 
fringe can be written as: 

P{Sf,(t)n,kn) = 

Ef=i PoHil + Ai sin[27rM/ + 0^J}, (3) 

where the summation is over all the N (target 
-I- field) stars within (and just outside) the FOV. 
In this quasi-monochromatic approximation, the 
fringes contributed by each field star will all have 
the same period, only their amplitudes and phases 



will be different. In that case, the final fringe will 
also be a pure sinusoid, and we can write it as 

P{Sf,4>n, kn) = 

P«{l + Ksin[27rA~(5/-|-0„]}, (4) 

where the parameters are related to those of the 
field star fringes by 

N 

P^ = E^oW' 

J = l 
N 

Pn Vn COS (j)n = ^ (n) A;^ COS 0^ , 

J = l 
N 

PnVnSm(j)n = ^ Pg^ (n) A{ siu </)^„ . (5) 

PnVn COS (pn and PnVnSiiKpn are respectively the 
cosine and sine components of the resultant fringe 
in each of the n narrow band channels. As as- 
sumed in Figure [5] (and without loss of general- 
ity), we take the target to be located at the center 
of the FOV, so that 9i = 0; its fringe phase (f)\ 
is therefore zero in all n channels, and the resul- 
tant phase (j)n is the confusion bias arising from 
the presence of the field stars. 

Note that Po{n), Am and are all functions of 
the nominal mean wavelength of the channel, the 
instrument parameters (as described earlier), and 
the intrinsic parameters of the stars such as their 
effective surface temperature T^f / , surface gravity 
log g, metallicity \og[Fe/ H], apparent magnitude 
rn;^, and location in the FOV. 

3.3.1. Recognizing confusion 

In each FOV we analyze, the model for SIM de- 
scribed above provides us with 80 measurements of 
the confusion bias If the brightness contrast 
between the target and the fields stars is large, 
we expect the channel-to-channel variation to be 
small and the resultant phase to be almost a lin- 
ear function of channel number, as described in 
i J2.2[ with the slope becoming zero if the target 
is located precisely at the origin of the adopted 
coordinate system in the FOV. For instance, the 
top panel in Figure [7] shows the phase spectrum 
for such a case; the target is a 10th mag A star, 
offset from the FOV coordinate origin by 1 mas, 
and is confused by a 15th mag M dwarf offset by 
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a (projected) angle of 100 mas. The linear part 
of the graph comes from the 1 mas offset of the 
A star, while the oscillatory part arises because of 
the presence of the faint M star. The latter be- 
comes relatively more important at longer wave- 
lengths, hence the amplitude of the phase "rip- 
ples" grows at smaller values of wavenumber. The 
middle panel shows another example of a binary, 
this time exhibiting a phase jump at short wave- 
lengths indicating a position offset of order 10 mas. 
The bottom panel of Figure [7] shows the phase 
spectrum for a crowded field; the fringe phase in 
this case varies erratically with channel number, 
but remains small. 

3.4. A confusion estimator 

Phase spectrum diagrams such as those just de- 
scribed (and the related delay spectrum described 
in H2.2p are useful indications of the possible pres- 
ence of a confusion bias in fields with relatively 
bright stars, but their varied character makes it 
less obvious how to combine them into a sin- 
gle measure which might be of value for fields of 
fainter stars. In a real observation with SIM, one 
answer to dealing with faint targets is to enlarge 
the bandwidth by binning channels on board the 
spacecraft before computing the fringe parameters 
and before transfer of the data to ground. On- 
board binning can be accomplished by adding the 
charge from several pixels as the CCD detector is 
read out. This has the advantage of reducing the 
relative contribution of read-out noise in the elec- 
tronics; the disadvantage is that information on 
the spectral variations of the fringes is no longer 
available, making it impossible to even recognize 
a confusion bias in a single observation^ A more 
thorough discussion of these issues is beyond the 
scope of this paper, but may be a useful topic for 
further work. 

In the framework of our simulations, we have 
developed an estimator for using all 80 channels 
of data based on summing the squared fringes in 
each channel after subtraction of the DC term. 
This quantity will vary as the fine delay Sf is 
scanned; minimizing this squared difference leads 
to a unique estimate of the mean astrometric de- 
lay, and hence the confusion bias. If the value of 



Averaging data in this way may also mask phase wraps 
even for unconfused targets. 




wavenumber"°(xl 0^) m 



Fig. 7. — Top: Phase spectrum to be expected in 
the case of a binary consisting of a mag 10 A star 
and a mag 15 M dwarf, separated by 100 mas along 
the direction of the baseline orientation. The re- 
sultant phase is a combination of a linear and an 
oscillatory function. Middle: A phase jump ap- 
pears when the binary is resolved at or close to one 
of the channels. Bottom: Phase spectrum to be 
expected in a crowded field. The resultant phase 
varies erratically with the wave-number. The mag- 
nitude of the phases and the resulting distribution 
across the channels determine the severity of the 
confusion bias. 




Fig. 8. — The confusion estimator as a function 
of fine-delay for the three cases shown in Figure [71 
The minimum value of and the corresponding 
fine-delay values are indicated for each case. There 
is some correlation between the amount of devia- 
tion of from zero to the amount of confusion 
bias; the exact nature of this relation is not inves- 
tigated here. 
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the estimator is zero at that delay to within the 
noise, the target is not significantly confused at 
that signal-to-noise level; non-zero values are an 
indication that the measurement is suffering from 
confusion. The estimator is written as 

n 

n 

where the delay spectrum is ^„ = 0„/27rA:„ (cf. 
W2.2\ and Figure [5]), and the other symbols have 
been defined previously in connection with equa- 
tion m As discussed in ^2.2[ ^„ is indepen- 
dent of the channel number n in the case of iso- 
lated targets (after adjustment for possible phase- 
wrapping). Thus, the quantity Sf+^n on the right 
side of Equation [H] can be considered as an effec- 
tive path- length delay of the target, and setting it 
to zero leads to the precise determination of the 
delay value S'j: at which becomes zero. This is 
therefore the delay offset of the target from the 
coarse delay position. Note also that at this value 
of delay, the fringes disappear, and the response 
pattern (Equation |4]) becomes the SED of the tar- 
get modulated by the throughput of the instru- 
ment. This discussion also shows that minimizing 
the estimator of equation [S] provides the cor- 
rect astrometric delay for an isolated target, and 
that the value of the estimator at that delay is 
exactly zero. 

From the discussion above it is also clear that, 
in the case of a crowded field, ^„ varies with chan- 
nel number and there is no single value of (5/ which 
will make each term of the sum in equation[6]equal 
to zero. The minimum value of will therefore be 
a nonzero positive number, so we will know there 
is confusion present. But the delay corresponding 
to that minimum is the estimator we will use for 
the target dclayQ In our simulations we know the 
target position, so we can calculate the true value 
of the target delay. We are therefore defining the 
confusion bias in the simulations to follow as the 
difference between the known delay of the target 
and the delay provided by minimizing the esti- 
mator. 



^To what extent this is an optimum estimator is presently 
uncertain. 



It should be emphasized that even though we 
have defined the confusion estimator in the con- 
text of these simulations, it would be possible to 
compute from real SIM data using the values 
of fringe parameters in each channel provided by 
the observations. A fictitious "fine delay" could be 
created and minimized as usual. To start with, 

could be estimated for targets known to be iso- 
lated and the noise floor could be determined. Any 
non-zero value of for other targets can then be 
used as an indication of confusion. 

4. Limiting values of confusion bias 

It is useful to have an estimate of the "worst 
case" values of confusion bias to be expected in 
a single SIM observation. While it is possible to 
imagine that the phasors for a specific distribu- 
tion of field stars could all add up "in phase" to 
produce a very large bias (cf. Figure [6|), this bias 
would be much reduced in the neighboring chan- 
nels; indeed, if there are many field stars, the ef- 
fects will be to simply raise the "noise level" of the 
average astrometric delay measurement on the tar- 
get. More serious will be those cases where only 
one or two relatively bright field stars are present. 
As an extreme example, we compute the confusion 
bias introduced by a single field star, as a function 
its (projected) angular distance from the target 
and the ratio of its brightness to that of the target. 
The SEDs of the field and target stars are assumed 
to be the same. At small separations, this bias will 
oscillate with a large amplitude and a period char- 
acteristic of the mean wavelength of the full pass- 
band. At larger separations the amplitude of the 
oscillations will fall owing to bandwidth apodiza- 
tion according to eauation lA6l Figure [5] shows the 
confusion bias (computed using the estimator 
as described in the previous section) as a function 
of projected separation for a field star with Am 
= -1-1. For SIM with 9.0 m (projected) baseline, 
a delay bias of 0.44 nm corresponds to an angu- 
lar position bias of 10 /las, which is approximately 
the expected single-measurement precision, so this 
value of delay bias is in some sense a "critical" 
value. This critical value is reached at a projected 
separation of 0.94" for a -1-1 mag field star. The 
critical projected separations for fainter field stars 
with Am = 2, 3, 4, and 5 mag are 0.17", 0.07", 
0.024", and 0.005" respectively. For instance if 
a field star is fainter by 3 magnitudes, it is not 
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0.0 0.5 1.0 1.5 

Projected separation 

Fig. 9. — Confusion-induced bias in the astrometric delay of a target star as a function of the projected 
separation of the 1 magnitude fainter field star, when both stars have identical SEDs. The inset indicates 
a magnified view of a portion of the plot. The dotted line in the inset indicates the delay bias of 0.439 
nm, which will correspond to a 10 /uas position uncertainty in an single observation. The critical projected 
distance beyond which the bias in the delay is less than 0.439 nm is 0.94 arcsec. 




02468 0246 

Magnitude difference (Am) Magnitude difference (Am) 



Fig. 10. — Confusion-induced bias in the astrometric delay of a target star as a function of the magnitude 
difference between the target and a fainter field star, both having identical SEDs. The biases for both close 
and wide binaries are plotted for a few arbitrarily-chosen projected separations. 
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likely to introduce a significant confusion bias as 
long as its projected separation from the target is 
more than 70 mas. Figure [10] shows the confusion 
bias as a function of magnitude difference for a few 
selected projected separations. 

5. Applications 

In this section, we present the results modeling 
confusion bias for a selection of fields representa- 
tive of several of the current set of Key Projects 
chosen for SIM subsequent to the first Announce- 
ment of Opportunity]^ As mentioned earlier, the 
important case of binaries as signal will be de- 
ferred to a later paper; here we consider such 
"companion" field stars as noise no matter where 
they are located in the FOV. 

5.1. Motions of nearby galzixies 

This project ("SIMDOG") aims at dynamical 
studies of the Local Group of galaxies by measur- 
ing their deviations from the Hubble fiowQ Ac- 
curate distances to the nearby galaxies from e.g. 
cepheids are a part of this work, but another es- 
sential component is the precise measurement of 
galaxy proper motions. Over the 5-year lifetime 
of the mission, repeated SIM observations of the 
brightest stars in galaxies out to distances of « 5 
Mpc are expected to yield these proper motions 
with sufficient precision. 

Initially, it might be assumed that proper mo- 
tion measurements of a galaxy with SIM would 
be insensitive to confusion by extraneous stars in 
the FOV. However, if the distribution of bright- 
ness over the FOV changes from one visit to the 
next, there is a possibility of confusion- induced bi- 
ases in the proper motions. One class of changes 
that has already been evaluated in the work of 
Raiagopal et al.l (|200lh involves biases caused by 
slight differences in pointing the SIM FOV on 
the target area during successive visits. A sec- 
ond class of changes in the brightness distribu- 
tion arises when Galactic foreground stars with 
their own proper motions are also present in the 
SIM FOV; this situation is likely to be avoidable 
merely by pointing SIM to a different target star 



'See the project descriptions at: 

http://planetquest.jpl. nasa.gov/SIM/sim^team.cfm. 
^The Principle Investigator for this project is E.J. Shaya. 



in the distant galaxy. The third class is not so 
easily circumvented; this arises from variations in 
the brightnesses of stars in the distant galaxy it- 
self. Variable stars are sufficiently common that 
one or more of sufficient relative brightness may 
occur within the FOV when SIM is pointed at a 
bright star in a nearby galaxy, and their variabil- 
ity may not be known ahead of time. Of course 
the target star itself may be variable, and this too 
may be initially unknown. It is this third class of 
biases we model here. 

5.1.1. Source model 

We assumed synthetic spectra of galactic su- 
pergiants of spectral type A and B with Te// = 
10000 K, log g = 2, and log[Fe/H] = -1.5. These 
are the most favorable stars for SIM observations 
in the V band, as described in the key project 
summary. 

5.1.2. Results 

The FOV is assumed to consist of the target 
star and one additional field star, both of which 
are stationary. The position of the target star is 
measured at two epochs, and the proper motion 
bias is estimated as the difference between the two 
positions. The target positions will contain some 
confusion bias owing to the field star and the gran- 
ularity of the background galaxy, but if everything 
remains the same the calculated proper motion 
will be zero. The brightnesses of the target and/or 
the field star are then changed, and changes in the 
confusion bias on the target position measurement 
will masquerade as a proper motion. 

We considered two cases: (a) the target and a 
field stars have the same brightness, but either one 
or both of them change in magnitude by 0.3 mag. 
between the two visits; and, (b) the field star is 
fainter than the target by 1 mag., but one of them 
changes by 0.3 mag. between visits. The bias in 
the proper motion for this case is a function of the 
magnitude difference and the angular separation 
(projected along the baseline) between the target 
and the field star. 

Our simulations indicate that there could be as 
much as « 800 /las of bias in the proper motions 
for galactic supergiants of spectral types A and B, 
in case (a) and « 300 /ias in case (b), if the field 
star is located within a fringe. If the projected sep- 
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aration of the field star is more than a fringe but 
less than 25 mas then the induced bias is « 30 /J,as 
in case (a) and « 20 iias in case (b). If the pro- 
jected separation of the field star is more than 50 
mas, then the bias in the proper motions owing to 
brightness variations of 0.3 mag is expected to be 
much less than the single measurement precision. 

5.1.3. Discussion 

A local group galaxy at a distance of 1 Mpc 
with a transverse velocity of 100 km/s would have 
a proper motion of « 20 /ias/yr, so that biases at 
the levels mentioned above could be problematic. 
Observing strategies to mitigate such biases will 
be necessary, such as repeating observations with 
changes of a few degrees in baseline orientation 
and/or observing several target stars in the same 
galaxy. As with the other cases of confusion bias 
which we have modeled here, the impossibility of 
obtaining a sufficiently-precise model of the bias 
precludes any simple correction for it; the best one 
can hope for is to be aware of the risks and plan 
the observations accordingly. 

5.2. Astrometric reference frame tie 

A list of K, 1300 Galactic stars spread more- 
or-less evenly over the sky will be repeatedly ob- 
served with SIM during the mission lifetime. The 
positions and proper motions of these "grid stars" 
will enable the definition of an astrometric refer- 
ence frame. Unfortunately, this frame may have a 
small residual rotation, a possible result of some 
non-random component to the measured proper 
motions of the constituent stars. This "roll" com- 
ponent will have to be removed by measuring the 
positions and (apparent) proper motions of 50 - 
100 distant quasars. These objects are generally 
too faint to be included in the regular astrometric 
grid program of SIM, but they can be measured 
to sufficient precision using longer integrations 

The quasars to be used in this study are dis- 
tributed all over the sky, and are relatively faint 
compared to grid stars. These quasars may there- 
fore be susceptible to confusion from faint (but 
numerous) A, K, and M stars in the Galactic fore- 
ground, especially near the Galactic Plane. It 
is this source of error in the astrometry of faint 



quasars which will be estimated here. The effects 
are exacerbated by the strongly differing SEDs of 
the quasars compared to the SEDs of faint red 
Galactic stars. 

5.2.1. Source model 

The models we adopt for this set of simulations 
consist of a target quasar, located for convenience 
at the center of the SIM FOV, and a field star lo- 
cated somewhere else within the FOV. Additional 
structures such as lumpy, extended jets, are not in- 
cluded. The quasar SEDs are chosen from a small 
but representative set of quas ar spectra from the 



•'First Bright Quasar Survey" (|White et al.ll2000l ) 



''The Principle Investigator for this project is K.J. Johnston. 



with redshifts ranging from to 3. These spec- 
tra have been extrapolated wherever they were 
not adequately defined over the wavelength range 
400 - 1000 nm by assuming a flat spectrum with 
values of the nearest known value; they are gen- 
erally brightest at the blue end of the spectrum. 
The Galactic star can be of spectral type AlV, 
KOV, or M6V. We do not include massive O and B 
type stars, since these stars are rare at the higher 
Galactic latitudes where the frame tie quasars are 
preferentially found. The foreground star is taken 
to be - 3 mag fainter than the target quasar, and 
located at PA = 5°, almost directly "above" the 
target. Four different radial separations between 
the quasar and the foreground star are modeled 
(25, 50, 100, and 200 mas), as weU as 18 different 
possible orientations (0° to 170° in steps of 10° ) 
of the interferometer baseline. 

5.2.2. Results 

The left panel of Figure [11] shows the single- 
measurement confusion bias as a function of the 
difference in brightness of the target and the loca- 
tion of the field star for a baseline position angle 
of 0°. The right panel of Figure [TT] shows a similar 
plot for a baseline position angle of 90°. Together, 
these figures show the typical bounds on the ex- 
pected biases. 

5.2.3. Discussion 

The simulations show a strong dependence on 
the SEDs of the quasar and the field star. For ex- 
ample, M stars cause significantly more bias then 
A stars of the same brightness, and there is a 
(smaller) dependence on the redshift of the quasar. 
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Fig. 11. — Absolute value of bias as a function of Am for field star separations of 25 (solid line), 50 (dotted 
line), 100 (dashed line), and 200 mas (dash-dot line) for AlV and M6V field stars. The PA of the field star 
is 5°. Left panel: Baseline PA = 0°. Right panel: Baseline PA = 90°. 



These biases can be reduced by choosing another 
baseline orientation as long as Am > 2 and the 
projected separation is > 50 mas for an A-type 
field star. The values for an M-type field star are 
3 mag and 100 mas, respectively. 

5.3. Taking the measure of the Milky Way 

SIM will make a major contribution to the 
study of the distributions of both dark and lumi- 
nous matter in the Galaxy through precision mea- 
surements of the distances and proper motions of 
different classes of Galactic starsj^ One of these 
studies involves measuring the distances to nearly 
100 bright M giants in the "Baade's Window" re- 
gion of the central bulge of the Milky Way. Unfor- 
tunately these are some of the richest star fields 
ever considered for SIM, and there are significant 
concerns about confusion-induced biases in single 
position measurements of these bulge stars. We 
have modeled typical examples of these fields and 
investigated whether judicious choices of the ori- 
entation of SIM's baseline can reduce these biases. 



The Principle Investigator for this project is S.R. Majewski. 



5.3.1. Source model 

We have constructed our models for these fields 
using IIST/WFPC2 imaging observations as a 
guide. For example. Figure [12] is an F555W ex- 
posure in the region of Baade's window (HST 
Proposal ID 8574; Target name FIELD180310- 
295143). This image has a scale of 0.0455" per 
pixel. We obtained photometric and astrometric 
information for more than 3800 stars in this im- 
age from J. Holtzman (private communication); 
his catalog included V & I magnitudes along with 
associated errors, X & Y position measurements, 
and the RA & DEC (J2000) of the stars. We ab- 
stracted a set of w 40 target FOVs from this im- 
age, each centered on a bright "target" star, and 
constructed an idealized model for each field con- 
sisting of a list of all stars lying within and just 
outside of the SIM FOV. The SED of each star was 
estimated as follows: First, we fitted a 2nd degree 
polynomial to a theoretical plot of V-I color index 
vs. effective temperature obtained from Kurucz's 
modelflll assuming that all the stars are solar-type 
dwarfs. This yielded a relation between the effec- 
tive surface temperature and the V-I color index 



^' \h,ttp://kurucz. harvard.edu/grids/gridPOOODFNEW/rijklpOOkOodfnew. dat\ 
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of the star. Using the known V-I colors for all 
stars in the catalog, we estimated the correspond- 
ing surface temperatures and hence the spectral 
types. A large fraction 85%) of the V-I col- 
ors lie in the theoretically-expected range of 0.5 - 
2.57; this is consistent with (but does not prove) 
the assumption that these stars have essentially 
zero metallicity. We also ignored any extinction 
corrections. All stars lying within a circle of 6" of 
the bright target star were considered as potential 
sources of confusion, since some light from these 
stars can be diffracted into the FOV. A baseline 
orientation was then chosen in the range 0° - 170° 
in steps of 10°, and the position of the target star 
calculated first without, and then including, the 
surrounding stars. The difference is the confusion- 
induced single-measurement astrometric bias. 

5.3.2. Results 

Of the 40 target fields we have simulated, 22 
show biases in excess of 10 /las if the same baseline 
orientation is chosen (90° for these simulations). 
The bias can be reduced below 4 //as in 39 out of 
40 targets simply by choosing a suitable baseline 
orientation. This suggests that some experimen- 
tation with the actual data set could therefore be 
useful, by taking the data at several baseline po- 
sition angles and rejecting anomalous points. Fig- 
ure[T3]shows an example of the astrometric bias as 
a function of baseline orientation for two different 
fields. For the first field (top panel), 12 out of 18 
baseline orientations show biases less than 10 /las. 
For the second field (bottom panel), 7 out of 18 
baseline orientations show correspondingly small 
biases. Also, note that the biases can occasionally 
be very large (e.g., at baseline PA = 20° in the 
second field). 

5.3.3. Discussion 

Our results generally show that, in these very 
crowded fields, the amount of bias grows with the 
number of field stars. This is expected, as con- 
sideration of the phasor diagram in Figure [6] will 
show. Furthermore, the largest contributions to 
the net bias are primarily caused by a small num- 
ber of judiciously-situated field stars. These stars 
can even be located outside the FOV. Figure [Ml is 
an example of such a case. In this field, there are 
56 field stars in our model. Two of them have pro- 
jected distances of only 2 mas and 13 mas, and are 




Fig. 12. — Image of the Milky Way Bulge 
(Baade's Window) retrieved from HST archival 
data. The image has a scale of 0.0455" per pixel. 
The image has been cropped and displayed as neg- 
ative. The black circle indicates the SIM FOV of 
3". 
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Fig. 13. — Single-measurement confusion-induced 
astrometric bias as a function of baseline orien- 
tation for models of two very crowded fields in 
Baade's Window. Note the logarithmic scale on 
the y axis; the biases can be large for some base- 
line orientations. The dashed horizontal line is 
drawn at a bias of 10 /ias, which is roughly the es- 
timated single-measurement error expected from 
noise and instrumental instability. 
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Fig. 14. — Left panel: One of the two field stars at position angle of 110° outside the field-stop, shown 
by filled dots, has a projected distance of about 2 mas and is fainter by 1.4 magnitudes. Middle panel: 
Excluding it from the model reduces the bias from 351 /ias to -2/ias. The other field star has a projected 
distance of about 13 mas, but it is fainter by 6.8 magnitudes; excluding it from the model makes little 
difference. Right panel: Somewhat surprisingly, a field star located well inside the FOV at a projected 
distance of 1.3" and brighter than the target by 0.95 mag also makes a negligible contribution to the net 
bias. 



fainter by 1.4 and 6.8 mag respectively, but they 
are both located outside the edge of the FOV and, 
at first sight, should not be a problem. Never- 
theless, the simulations show that the measured 
target position is biased by 351 ^as. When the 
field star at 2 mas is excluded, the bias drops from 
351 fias to -2 /ias. Further excluding the 13 mas 
field star makes very little difference; this star is 
probably too faint to cause problems, or it is at a 
fortuitous location on the fringe pattern. There is 
a third field star which would seem a likely candi- 
date for causing a bias; it is brighter by 0.95 mag 
and located inside the FOV, at a projected dis- 
tance of 1.3". However, excluding this star makes 
little difference to the bias, as is shown in the right 
panel of Figure [TH 

5.4. Cluster distances 

This project is aimed at measuring the dis- 
tances to a set of selected open and globular clus- 
ters in the Galaxy in order to determine their pre- 
cise ages by combining SIM parallax data with 
other existing datal^ A few stars from each of 
these open and globular clusters will be chosen as 
targets for SIM astrometry. However, these tar- 
gets are necessarily in crowded stellar fields, and 
the examples of the previous section suggest that 
the confusion-induced biases could be large. We 
have modeled these biases in the typical case of 



The Principle Investigator for this project is G. Worthey. 



NGC 6440, an old globular cluster at a distance of 
« 8.4 kpc, and located near the Galactic center at 
1 = 7.7°, b = +5.8°. NGC 6440 is a well-studied 
globular cluster, with a core radius of 8". It is 
one of the 14 globular clusters known to have X- 
ray emission. Specifically, we have addressed the 
following questions with our models: 

• What is the typical confusion bias for the 50 
brightest targets in NGC 6440? 

• At what area density of bright field stars are 
the target biases tolerable? 

5.4-1. Source model 

Our model for NGC 6440 is based on an HST 
image obtained with the Planetary Camera of 
WFPC2, as shown in Figure [15] and kindly pro- 
vided by G. Worthey. Most of the bright stars 
in this picture are K giants; we generated a cat- 
alog of these stars using the 'SExtractor' soft- 
ware Although it is well known that this ap- 
proach does not provide accurate photometry in 
crowded fields, our goal here is not to obtain pre- 
cise corrections to the astrometric observations 
(which would require micro arc-second positions 
as well as accurate magnitudes) but rather to elu- 
cidate the general level and nature of the biases 
which may be expected. The number of sources 



' Bertin, E., 

\http://terapix.iap.fr/rubrique.php?id.rubrique=91\ 



17 



Fig. 16. — Number density of stars as a func- 
tion of radial distance from the core of the cluster 
NGC 6440. The magnitudes of the stars lie in the 
range from 16.19 to 25.52 



Fig. 15. — PC image of the Globular Cluster 
NGC 6440 cropped and displayed as a negative. 
The image scale is 0.046 "/pixel. The black circle 
indicates the SIM FOV of 3". 



sumed that all the stars have solar metallicity. The 
spectra so determined were renormalized to ac- 
count for the distance to the cluster and the mag- 
nitude of the field star. 



which can be obtained with this software depends 
on the assumed value of the background intensity; 
after a few trials, we settled on a background value 
of 1.5 X arms, which yielded about 2400 stars with 
positions and V magnitudes. A visual compari- 
son of the selected objects with the original image 
showed a good correspondence. The area density 
of these stars as a function of the radial distance 
from the center of the cluster is shown in Figure 
[TCI We then selected 51 circular fields at different 
locations in the cluster (one is shown in Figure fTS)) : 
each field is centered on a bright star (presumed 
to be a K-giant). The model constructed for that 
field consists of a list of the positions and magni- 
tudes of the target star and all field stars within a 
radius of 3" of the target; this allows for diffraction 
of hght into the 3" (diameter) SIM FOV from stars 
which are located beyond the edge of the aperture. 

The SEDs of the model stars were determined 
from estimates of surface temperature, gravity, 
and metallicity, as follows. From the V mag- 
nitudes (with the assumption that all stars are 
K giants) we attributed a temperature to a star 
randomly in the range 4000-4750 K (4000, 4250, 
4500,4750) and a surface gravity randomly in the 
range of 1.5 to 4 with a step size of 0.5. We as- 



5.4-2. Results 
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Fig. 17. — Probability that the confusion bias ex- 
ceeds 10 /xas in a single measurement as a function 
of radius of the globular cluster NGC6440. 

In each of the 51 fields we obtained the con- 
fusion bias for 18 different orientations (position 
angles) of the SIM baseline (0° to 170° in steps of 
10°). The fields used in this set of simulations typ- 
ically have target stars at V « 16-18, with bright- 
est field stars at V « 16 or fainter, down to « 25 
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(with some exceptional cases; in 5 out of 51 fields, 
brightest field stars were brighter than the target 
by 1.2, 1.6, 3.2, 2, 4.6 magnitudes; 3 of them were 
within the core radius and the remaining two were 
within 16 arcsec radius.). Of aU the cases we ana- 
lyzed in each field, 20% of them showed biases in 
excess of 10 /ias in all 18 baseline orientations con- 
sidered; Figure [17] shows the probability of a sin- 
gle measurement confusion bias exceeding 10 fxas 
as a function of radius of the cluster. For a given 
field, if the confusion bias was more than 10 fj,as 
in 12 out of 18 orientations, the probability was 
estimated as 2/3 at the radial distance of the field 
from the center of the globular cluster. This fig- 
ure indicates that, for NGC6440, the probabihty 
of having large confusion bias can be as large as 
95% even at the radius of 15" and there could 
be as much as 40% probability of large confusion 
bias at large radial distances. The probability of 
obtaining large confusion bias drops at larger ra- 
dial distances in the cluster, as the local stellar 
surface density decreases. This latter point is il- 
lustrated in Figure [m which shows the confusion 
bias as a function of baseline orientation for two 
different fields, located at radial distances of 3.6" 
and 20" from cluster center. For the field at 3.6" 
(which is actually inside the cluster core), the bi- 
ases are always large (note the logarithmic scale) 
at all baseline orientations. For the field at 20", 
the error is larger than 10 /ias in 30% of the base- 
line orientations. 

5.4-3. Discussion 

The statistics of these results indicate that 
the probability of having an unacceptable level 
(> 10 /ias) of confusion bias in the positions of 
target stars is will exceed 30% for a single obser- 
vation carried out at any random baseline orien- 
tation if the stellar density in the area exceeds 0.4 
stars per square arc-second (compare Figures [12] 
and I17[ In other words, about 3 stars within the 
3" FOV is acceptable). 

However, the focus of this particular Key 
Project is on distances to objects, and hence on 
the parallaxes of the target stars. How do our 
results apply to the measurement of parallaxes? 
Such a measurement could be carried out with 
SIM in a special way, namely, ensuring that the 
baseline orientation was identical for observations 
taken « 6 months apart in time. This is per- 
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Fig. 18. — Confusion bias as a function of base- 
line orientation for model fields located at 3.6" 
(top panel) and 20" (bottom panel) from the cen- 
ter of NGC 6440. Target fields at the periphery of 
clusters are clearly preferable, but may still have 
an appreciable probability (here ~ 30%) for a sub- 
stantial confusion bias. 

haps possible in principle, but very unlikely in 
practice. First, the baseline orientation can not 
easily be set to « 1°; indeed, it is in fact gen- 
erally not even known until the baseline vector 
has been calibrated. This means that parallaxes 
will suffer from the same magnitude of confusion 
bias as do regular position measurements. In the 
case of the cluster Key Project, a possible strat- 
egy could be to choose a number of target stars 
each at some large distance from the cluster cen- 
ter (such that the stellar density is lower than the 
limit cited above), and obtain parallaxes on each 
of them using the standard SIM observing strat- 
egvF^ Discordant single-measurement data points 
would then simply be rejected from the data set 
as a likely consequence of confusion. 

6. Limitations of the models 

The first, and perhaps most painful, limitation 
of the modeling we have described here is that, in 
spite of our apparent ability to compute the con- 
fusion bias in any single measurement with high 
precision, the results are in general not likely to be 
sufficiently accurate to provide actual corrections 
to the target positions. The simple reason for this 



This strategy is still under development and includes a com- 
plex calibration program. 
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is that the true field star positions are not known 
with accuracies of the same order as the expected 
single- measurement precision of SIM, w 10 /xas. 
Whether the modeling might be sufficient in any 
particular case will depend on the specific details 
of that field on the sky; our modeling tools (de- 
scribed in the next section) can then be used to 
evaluate the biases by e.g. varying the positions 
and brightnesses of the field stars over plausible 
ranges. 

The second limitation of our modeling concerns 
the instrument model we have adopted for SIM. 
There are several simplifications we have made 
which could be problematic: 

1. In our numerical simulations, we have as- 
sumed that the FOV is exactly circular, and 
that all the photons diffracted into it will be 
collected (with some efficiency) by the de- 
tector. However, in practice the focal plane 
of SIM's camera will be covered by a CCD 
detector with pixels of « 2" on a side, and it 
is presently planned to average the data in 
three neighboring rows. 

2. We have assumed that the photons collected 
in each FOV will be dispersed in the camera 
in some way into 80 channels with central 
wavelengths and bandwidths as specified in 
Appendix lB.il In fact, a thin prism will be 
inserted into the light path before imaging 
onto the CCD detector, turning the camera 
into an objective prism spectrograph. Such 
instruments require special calibration and 
are subject to a degree of internal confusion 
caused by overlap of spectra from field stars 
with the spectrum of the target star. 

3. We have estimated the fringe parameters 
(total power, fringe amplitude and fringe 
phase) analytically, and neglected the details 
of just how they will be measured on-board 
the spacecraft. 

4. We have assumed that the throughput of the 
system is the same for different locations in 
the FOV. In fact, the presence of the prism 
will cause the throughput (and the disper- 
sion) to change depending upon the angle of 
incidence at the prism, and hence to be dif- 
ferent for the target and for the field stars. 



We have carried out a further study using a 
more sophisticated instrument model which re- 
moves the first three of these limitations. The 
details of this extended model will be described 
elsewhere ( Sridharan Sz Allen l2007d) . but we give 
here the results of repeating the bias computations 
on a subset of the fields we have described in the 
present paper. Table [T] compares the biases ob- 
tained using the "simplified" approach presented 
in this paper with those obtained using the more 
detailed model of the focal plane. The cases listed 
are taken from the Quasar Frame Tie key project 
(cases 1-5) and a selection of binary models (cases 
6-8). Perhaps not surprisingly, the differences 
are roughly proportional to the magnitude of the 
computed bias, although these differences are gen- 
erally at a level of only a few percent. Unfortu- 
nately, the utility of the results from this more so- 
phisticated instrument model is still compromised 
by inaccurate input data on the field stars, and 
since it also comes with considerable additional 
complexity we have not used it further in this pa- 
per. 

7. Modeling Tools 

The source and instrument models described 
here have been implemented in a suite of programs 
written in the IDL programming system. These 
programs are availablqlj for further experimenta- 
tion, including the code developed for estimating 
confusion bias in the specific target fields described 
in this paper. 

8. Summary and Discussions 

We have examined the bias that can occur in a 
single measurement with SIM owing to the pres- 
ence of field stars within the FOV. In order to 
accomplish this task we have presented a model 
for the SIM interferometer, and a description of 
how SIM carries out a single measurement of the 
position of an isolated target star. The measured 
instrument response is then perturbed by adding 
a field star to the model FOV; the difference in 
the angles measured in the two cases is called the 
"confusion bias" . The extremes of this bias are 
calculated for the specific (but common) case of a 
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Table 1 

Comparison of the single-measurement confusion bias obtained with the present simplified 

MODEL AND WITH THE MORE DETAILED MODEL OF ISRIDHARAN fc ALLEnI (|2007bI '). 



Case 


Model 


Simple 


Detailed 


Difference 


Number 


parameters 


(/I as) 


(^as) 


(/las) 


1 


Q(z=0), AlV, 2, 50, 5, 90 


442 


431 


11 


2 


Q(z=0), M6V, 2, 50, 5, 90 


1072 


1036 


36 


3 


Q{z^2), AlV, 2, 50, 5, 90 


224 


215 


9 


4 


Q(z=2), M6V, 2, 50, 5, 90 


614 


592 


22 


5 


Q(z=2), M6V, 2, 50, 5, 10 


-46 


-46 





6 


AlV, BIV, 3, 1500, 10, 90 


-1 





1 


7 


AlV, M6V, 2, 50, 90, 90 


-59 


-61 


2 


8 


AlV, M6V, 2, 25, 90, 90 


-53 


-46 


7 



Note. — In column 2, Q(z=:0), AlV, 2, 50, 5, 90 means target quasar with 
redshift 2, AlV field star, Am — 2, radial distance 50 mas, PA = 5°, baseline 
orientation 90°. 



binary system in order to illustrate its main prop- 
erties. 

A number of source models are then developed 
which resemble the fields to be studied by SIM 
in several of the Key Projects already selected for 
inclusion in the initial mission science program. 
An unconfused version of the source model con- 
sisting only of the main target star is used as a 
reference measurement, and the results compared 
with a measurement made on the fully-populated 
field. The difference is the confusion bias in a sin- 
gle SIM measurement. Observations are simulated 
at various orientations of the interferometer base- 
line, and variants of the full field are examined in 
order to understand the sensitivity of the bias to 
structural details in the field. 

The magnitude of the confusion bias is found 
to depend on a number of factors, some obvious, 
others perhaps less so: 

• the relative brightnesses of the target and 
the field stars; 

• the shapes of the SEDs of the target and field 
stars; 

• the angular separation of the stars from the 
center of the FOV; 



• the angular separation of the field stars from 
the target star as projected on the interfer- 
ometer baseline; and, 

• the baseline orientation. 

The largest contributions to the confusion bias in 
a crowded field come from a small number of stars 
having small projected angular separations from 
the target, but these stars may actually be located 
outside of the FOV. Field stars which are less than 
4 mag fainter than the target and which have pro- 
jected separations within 100 mas of the target are 
potentially the most troublesome. 

The results of this study provides the under- 
standing and the tools required to examine the 
likelihood of confusion bias in any single measure- 
ment with SIM. Unfortunately, data on the field 
stars in any specific FOV (especially their posi- 
tions) is not likely to be available with sufficient 
accuracy to actually remove this biasj^ Our study 
nevertheless suggests some strategies for recogniz- 
ing the presence of confusion bias and for deal- 
ing with it, both in the observation planning stage 
and in the data reduction stage. These strategies 



There are some possible exceptions one could imagine, but 
we have not explored them further here. 
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might include the following: 

• While dealing with crowded fields, avoid 
fields with star densities in excess of 0.4 stars 
per square arcsec. 

• If avoidance is impossible, evaluate the like- 
lihood of confusion in the field by using the 
tools developed here. 

• If confusion is likely, try to reduce your sen- 
sitivity to it by planning the observing pro- 
gram so that data is taken at the least sensi- 
tive orientations of the interferometer base- 
line. 

• If too little is known about the specific field, 
plan to distribute the available observing 
time over a number of orientations of the 
interferometer baseline which differ by a few 
degrees from each other. Inconsistent val- 
ues in the data set can then be rejected with 
motivation. 

• If confusion is suspected in a given set of 
observations for which no prior data exists, 
acquiring new imagery from e.g., speckle or 
adaptive optics imaging would be useful for 
building a model. 

There is one additional strategy suggested by our 
more accurate model of the SIM focal plane. The 
CCD detector in the focal plane of SIM's camera is 
planned to have pixels which are smaller than the 
diameter of the FOV. If the data in the individual 
pixels can be made available, it would be possi- 
ble to choose e.g. only the central pixel, thereby 
effectively reducing the FOV and possibly attenu- 
ating an offending field star. The penalty of fewer 
target photons could then be offset by a reduc- 
tion in the level of confusion bias. This possibility 
will be discussed in more detail in a future paper 
(ISridharan fc Allenll2007bh . 

We wish to emphasize that the results of this 
paper refer to a bias which may be present in a 
single measurement of angular position with SIM. 
The determination of the full set of astrometric 
parameters (position, parallax, proper motion) on 
any SIM target will be done with a number of 
measurements, reducing the effects of any single 
anomalous point. Furthermore, the ultimate ac- 
curacy of the results depends on an extensive cal- 
ibration program to determine the instrumental 



parameters, including the baseline length and ori- 
entations for each field observed. 
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A. Michelson interferometer response 

Michelson interferometers are used in astronomy at wavelengths from the radio to the optical for spec- 
troscopy, for astrometry, and for synthetic imaging. The major equations giving the response of such inter- 
ferometers for the latter two applications are summarized here. 



A.l. Interferometers for astrometry and imaging 

At any given wavelength A and baseline separation Bij (both e.g. in meters), the response of a simple 
adding Michelson interferometer of the ty pe used for as t romet ry and imaging of partially coherent light in 
optical astronomy can be modeled as (e.g. Born fc WoH ( 1975[ ). p. 491 et. seq.): 



P(Ay,A:) = Po[l + ^cos(2^A:A,,)], (Al) 

where Pq is proportional to the target brightness (e.g. in photons/sec), A < 1 is called the fringe amplitude 
and depends on the target structure, k = 1/A is the wavenumber (in e.g. meters"^), and Ay is the total 
optical path difference OPD (e.g. in meters) between the two sides (or "arms" ) of the interferometer. For point 
source targets, A — 1, and the normalized version of equation lAll can be considered as the interferometer's 
far-field "point spread function" (PSF), similar in concept to the PSF of a fiUed-aperture telescope. 

The total OPD Ay consists of components on each side i and j of the interferometer, and each of those 
consists of an internal and an external part. Furthermore, some of the internal delay is "fixed", and some 
may be "variable". Referring to Figure [TOl we can write: 

Ay — [total side j delay] — [total side i delay] , 

where the components on each side are: 



side j internal delay = fixed delay -I- variable delay 

= dj -\- (5, and 

side j external delay — by construction. Further, 

side i internal delay — fixed delay -|- variable delay 

— -|- by construction, and 

side i external delay — Bij sm9. 

This last equation is simple geometry, as Figure [19] shows, but it is often stated as the dot product of a vector 
B parallel to the direction of the baseline with a (unit) vector D in the direction for which we wish to compute 
Ay . D makes an angle of a with B, see Figure [TH Then the "side i external delay" can be written as B • 
D — Bij cos a. Defining a = 9 + tt/2 we have "side i external delay" = Bij sinfl, with 9 defined as shown 
in Figure [191 in fact, 9 is the half-opening angle of a cone with axis parallel to the interferometer baseline. 
Note that 9 is also the angular distance of the direction of interest from a direction perpendicular to the 
baseline (— 7r/2 < 9 < -|-7r/2). Our idealized interferometer has a constant response in directions orthogonal 
to the baseline; the actual field of view will be further restricted by the practicalities of its design. Note 
further that equation I A 1 1 also closely describes the pattern in transmission, such as would occur if a laser (or 
CW radio transmitter) would be sent from an appropriate point in the beam path backwards through the 
interferometer. This is an expression of a general reciprocity theorem, familiar in radio engineering, which is 
itself rooted in the symmetry of the wave equation for the electromagnetic field with respect to the direction 
of time. 

Finally, equation lAll applies only at a single wavelength; one must also account for the finite bandpass. 
In essence the response becomes a sum of many different patterns within the bandpass, and the coefficient of 
the "cos" term in equation [XT] will be multiplied by a factor which is essentially the Fourier Transform of the 
band shape for that specific channel. Without detailed knowledge of that band shape we can not calculate 
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Fig. 19. — Sketch of the basic Michelson interferometer as used for astrometry and imaging in radio and 
optical astronomy. B is a vector in the direction of the basehne, and D a unit vector in the direction in 
which we want to compute the interferometer's response. 6 = a — Tr/2, and Pq, A, and are defined in the 
text. 
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the final form of equation lAll however, the results are likely to be amenable to a simple parameterization, 
using the following definitions. The mean wavenumber kn for channel n is defined as: 

[kK{k)dk 

jF^{k)dk ' ^'^'> 

where Fn{k) is the band shape of channel n expressed in wave numbers {k = 1/X = v/c, where v is the 
frequency in Hz). The coherence length A„ (e.g. in meters) for channel n is defined here as: 

A„ « 1/(AA:)„ (A3) 

where (Afc)„ is the full width at half maximum of an assumed Gaussian band shape Fn{k) centered at fc„. 
The coherence length is the scale size of the wave packet formed by the collection of photons at various 
neighboring wave numbers that make up the broad-band signal. It is given by the expression A„ = A /AA. 
With these definitions we can write our model expression for the fringe pattern in channel n as: 

PiA^jXiAn) = Pq{1 + Ae'^i^) cos(2^i;:A,y)} (A4) 

where C = 7r^/41n2 = 3.56. Suppose further that we design the interferometer such that d,; — dj (or define 
6 = such that this occurs), then the OPD becomes: 

Aij =6- Bij sine. (A5) 

Inserting this last relation into equation IA4I gives the final (approximate) expression for the interferometer 
response in a single channel of finite bandwidth: 

r S-B^j sin 9 1 2 

P{S,e,k^) = Pq{1 + Ae^^-^^i ^ i cos(27rA~[(5-B,jSin6i])}. (A6) 

Considered as a function of 9 this equation describes the fringe pattern anywhere in the sky for a given value 
of internal delay S. Alternatively, equation I A6I describes the response to a source located at a fixed 6 in the 
sky as a function of the internal delay S. 

Equation lA6l needs to be modified in order to render it specific to SIM. First, SIM's optical beam combiner 
effectively adds an additional 7r/2 of phase delay to one of the incoming beams, turning the "cosine" function 
into a "sine" . Second, SIM always "observes" in a direction closely perpendicular to the baseline orientation 
(cf. Figure fT9|) . so that 6 is small and we can set sin 6* w 9. Third, the exponential term describing the 
amplitude decorrelation can be set to unity, as the bandwidths of the individual channels are narrow. With 
these changes, and ignoring the subscripts ij of the baseline _B, we can write equation IA6I for any channel 
with mean wavenumber k as: 

PiS,e,k) = Po{l + A8m{27rk[S - B9])}. (A7) 



B. Instrument and Model Details 

B.l. Current SIM parameters 

Our SIM instrument model assumes the following values for the various parameters required in the sim- 
ulation of the interferometer response: 

Baseline length: 9.000 m for the astrometric interferometer. 

Collector size: Siderostats, outer diameter 304.5 mm, inner diameter 178 mm, net area = 479.378 cm^. 
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Field stop size: Nominally 3" diameter. However, this can be chosen in the simulation code to be any one 
of four possible values (1,2,3 or 4"). 

Number of channels and bandwidth: The design for the fringe disperser in SIM has 80 narrow-band channels 
with bandwidths progressively increasing from 1.8 nm at 401.9 nm, to 24.9 nm at 985.5 nm. The band 
shape is likely to be approximately Gaussian, and this is what we have used; however, the simulation 
code allows the user to choose a rectangular passband instead. The list of central wavelengths and 
bandwidths assumed for these simulations can be found at the author's web siteF^ 

Throughput (Th): The Throughput varies across the 80 channels, depending on the reflection/transmission 
optics and the QE and spectral response of the detector. We found that the experimentally-measured 
values could be modeled best using the relation Th = Y1\=q c^A' with co — 18.2, ci — —338.5, C2 = 
2076.2, C3 = -6239.4, C4 = 10463.3, C5 = -10029.1, cq = 5151.4 and C7 = -1102.1, for A expressed in 
microns. 

Pointing accuracy: There are two parts to the "pointing" of SIM: The angle tracker will center the target 
in the instrument FOV and superpose the images from each side of the interferometer to within 10 mas 
for bright targets Vtarget < 15, and 30 mas for faint targets 15 < Vtarget < 19. The fringe tracker will 
set the coarse delay on the target with an accuracy of 10 nm (for Vtarget < 10, corresponding to a 
maximum offset of « 0.2 mas at 500 nm. 

Spectral dispersion: A thin prism is inserted into the light path after beam combination, turning the 
instrument into an objective prism spectrograph, such that the images of stars in the focal plane are 
stretched out into spectra in a direction approximately parallel to the projection of the SIM baseline 
on the sky. For a target at the center of the FOV, the experimentally-measured dispersion can be 
modeled with the polynomial x = X]^=o '^i^"' with cq — 5.1, ci = —2.7, C2 = 1.0, C3 — —0.23, C4 — 0.03, 
C5 = 0.002, and cg = 0.00005, for v in units of 10^'^Hz and x in mm. 

Focal plane camera: The camera has pixels of size 24 fi aligned along the direction of dispersion; however, 
we ignore the pixellation of the focal plane camera in the simulations described in this paper. 

Overall measurement precision: The design goal is to make a single measurement of the angular position of 
a target on the sky with a precision of < 10 fias. As the fringe period X/B is typically 10 - 20 mas and 
the baseline length is about 10 — 20 x lO^A, the design requirement on the single-measurement angular 
precision corresponds to w 1/1000 of a fringe. 

B.2. Estimation of Total Power 

In this Appendix we present the mathematical details of how we estimated the total light from the target 
and field stars that lie within and just outside the 3" FOV of SIM. 

For the target star, the total energy received by the interferometer in each channel is multiplied by the 
throughput and the encircled energy within the aperture at the mean wavelength A (to account for the FOV 
set by the field-stop) and integrated over the bandpass to obtain the total power in that channel. It is then 
multiplied by hc/X to obtain the number of photons/second/channel. 

In the case of the field stars, the total energy is obtained as follows: Consider a Cartesian coordinate 
system centered on the target star which is also at the center of the FOV. Suppose that there is a field 
star at (xo,j/o)- The energy of that part of the diffracted image that lies within the field stop is obtained 
by calculating the integral of the product of the point spread function (PSF) centered at {xo,yo) and a 
pill-box function representing the field-stop. This is equivalent to integrating the product of the on-axis PSF 

^' \http://www.stsci.edu/^rjallen/sim/\ 
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Fig. 20. — Fraction of the energy recorded for a star located at a given radial distance from the center of 
the (circular) FOV of SIM, for field stops of 2, 3 and 4" diameter, and at several different wavelengths. 

representing the field star with a shifted version of the pill-box function centered at (xq, yo)- The expression 
is: 



and EF is the estimated flux, PSF is the point spread function, and FSD is the field stop diameter. From this 
representation it is clear that the energy of the field star that lies within the field stop is the convolution (or 
correlation) of the on-axis PSF of the field star with a pill-box function located at (a;o,yo)- It is convenient 
to normalize this correlation function and obtain a fractional number which, when multiplied by the total 
energy received by the interferometer from the field star, provides the desired value. This value is further 
multiplied by the throughput of the system and converted into photons/second/channel for each field star. 
Figure [20] shows this overall transmission function for several possible field stop sizes and wavelengths. We 
have computed this response function at the mean wavelength of each channel. 




(Bl) 



where the "Circ" function is defined as: 




27 



REFERENCES 

Bodcn, A. 2004, "Space Intcrfcromctry Mission, 
Data Product Document," Jet Propulsion Lab- 
oratory & Michelson Science Center, Caltech, 
Pasadena, California. 

Born, M., & Wolf, E. 1975, "Principles of Optics," 
(Oxford, Pergamon) 

Dalai, N., & Griest, K. 2001, ApJ, 561, 481 

Gaskill, J.D. 1978, "Linear Systems, Fourier 
Transforms, and Optics" (New York, Wiley), 
Section 2-2 

Hecht, E. 2002, "Optics" (Addison Wisely) 

Milman, M. H., & Basinger, S. 2002, Appl. Opt., 
41, 2655. 

Milman, M. H., & Truyshev, S. V. 2003, Opt. 
Eng., 42, 1873. 

Phillips, E.G. 1961, "Functions of a Complex Vari- 
able" 8tli ed. (Edinburgh, Oliver k Boyd), 5 

Rajagopal, J., Boker, T., & Allen, R.J. 2001, 

PASP, 113, 1232 

Sridharan, R. & Allen, R. J., 2007b, in preparation 

Takeuchi, T., Velusamy, T., & Lin, D. N. C. 2005, 
ApJ, 618, 987 

White, R.L., Becker, R.H., Gregg, M.D., Laurent, 
M.S.A (& and 11 co-authors) 2000, ApJS, 126, 
133 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v5.2. 



28 



